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ABSTRACT 

We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet 
ft Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve 
high signal-to-noise, high precision measurements. While /J Pic b’s near-IR colors mimic those of a standard, 
cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few 
substellar/planet-mass objects — i.e., k And b and 1RXJ 1609B — match ( J > Pic b’s JHK S L' photometry and its 3. 1 /i m 
and 5 jtm photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and 
large (~60 /im) dust grains fail to reproduce the /I Pic b spectrum. However, models incorporating thick clouds 
similar to those found for HR 8799 bcde, but also with small (a few microns) modal particle sizes, yield fits consistent 
with the data within the uncertainties. Assuming solar abundance models, thick clouds, and small dust particles 
((a) = 4 fim), we derive atmosphere parameters of log(g) = 3.8 ± 0.2 and 7' c ir = 1575-1650 K, an inferred mass 
of 7 44 My, and a luminosity of log(L/L 0 ) ~ —3.80 ± 0.02. The best-estimated planet radius, «G.65 ± 0.06 R j, 
is near the upper end of allowable planet radii for hot-start models given the host star’s age and likely reflects 
challenges constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with 
hot-start model predictions if /3 Pic b is younger than ~7 Myr, consistent with a late formation well after its host 
star’s birth ~12 4 4 Myr ago. 
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1. INTRODUCTION 

The method of detecting extrasolar planets by direct imaging, 

even in its current early stage, fills an important gap in our knowl- 

edge of the diversity of planetary systems around nearby stars. 
Direct imaging searches with the best conventional adaptive op- 
tics (AO) systems (e.g., Keck/NIRC2, Very Large Telescope 

( VLT) /NAOS-CONICA (NaCo), Subaru/HiCIAO) are sensi- 
tive to very massive planets (M >5-10 My) at wide separation 

( a ~ 10-30 AU to 100 AU) and young ages (t < 100 Myr), which 

are not detectable by the radial velocity (RV) and transit methods 

(e.g., Lafreniere et al. 2007a; Vigan et al. 2012; Rameau et al. 

2013; R. Galicher et al. 2013, in preparation). Planets with these 
masses and orbital separations pose a stiff challenge to planet 
formation theories (e.g., Kratter et al. 2010; Rafikov 2011). 
Young, self-luminous, directly imageable planets provide a crit- 
ical probe of planet atmospheric evolution (Fortney et al. 2008; 
Currie et al. 2011a; Spiegel & Burrows 2012; Konopacky et al. 
2013). 

The directly imaged planet around the nearby star /> Pictoris 
(ft Pic b) is a particularly clear, crucial test for understanding 
the formation and atmospheric evolution of gas giant planets 
(Lagrange et al. 2009, 2010). At 12 4 4 Myr old (Zuckerman 
et al. 2001), the ft Pic system provides a way to probe planet 


atmospheric properties only «s5— 10 Myr after the disks from 
which planets form dissipate (»G-10 Myr; e.g., Pascucci et al. 
2006; Currie et al. 2009). Similar to the case for the HR 8799 
planets (Marois et al. 2010b; Fabrycky & Murray-Clay 2010; 
Currie et al. 2011a; Sudol & Haghighipour 2012), ft Pic b’s 
mass can be constrained without depending on highly uncertain 
planet cooling models: in this case, RV-derived dynamical mass 
upper limits, when coupled with the range of plausible orbits 
( a ~ 8-10 AU), imply masses less than ~ 1 0— 1 5 My (Lagrange 
et al. 2012b; Currie et al. 201 lb; Chauvin et al. 2012; Bonnefoy 
et al. 2013), a mass range consistent with estimates derived from 
the planet’s interaction with the secondary disk (Lagrange et al. 
2009; Dawson et al. 2011). 

Furthermore, while other likely/candidate planets such as 
Fomalhaut b and LkCa 15 b are probably made detectable by 
circumplanetary emission in some poorly constrained geometry 
(Currie et al. 2012a; Kraus & Ireland 2012), ft Pic b’s emission 
appears to be consistent with that from a self-luminous planet’s 
atmosphere (Currie et al. 2011b; Bonnefoy et al. 2013). Other 
objects of comparable mass appear to have formed more like 
low-mass binary companions. Thus, combined with the planets 
HR 8799 bcde, ft Pic b provides a crucial reference point with 
which to interpret the properties of many soon-to-be imaged 
planets with upcoming extreme AO systems like GPI, SCExAO, 
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Table 1 

Observing Log 


Currie et al. 


UT Date 

Telescope /Instrument 

Mode 

Pixel Scale 
(mas pixel -1 ) 

Filter 

tint 

(s) 

N images 

AP.A. 

(deg) 

2011 Dec 16 

VLT/NaCo 

Direct 

13.22 

J 

50 

48 

26 

2012 Jan 1 1 

VLT/NaCo 

Direct 

13.22 

H 

40 

92 

34.6 

2012 Dec 15 

VLT/NaCo 

Direct 

27.1 

M' 

20 

176 

70.2 

2012 Dec 16 

VLT/NaCo 

Direct 

27.1 

l! 

30 

112 

67.7 

2012 Dec 16 

VLT/NaCo 

Direct 

27.1 

[4.05] 

30 

64 

57.4 

2012 Dec 23 

Gemini/NICI 

Direct 

17.97 

[3.09] 

38 

60 

30.2 

2012 Dec 26 

Gemini/NICI 

Direct 

17.97 

[3.09] 

38 

60 

31.7 

2013 Jan 9 

Gemini/NICI 

O'.' 22 mask 

17.97/17.94 

H/K s 

11.4 

117 

41.1 


and SPHERE (Macintosh et al. 2008; Martinache et al. 2009; 
Beuzit et al. 2008). 

However, investigations of /I Pic b’s atmosphere are still 
in an early stage compared with those for the atmospheres 
of the HR 8799 planets and other very low-mass, young 
substellar objects (e.g., Currie et al. 201 la; Skemer et al. 2011; 
Konopacky et al. 2013; Bailey et al. 2013). Of the current 
published photometry, only K s (2.18 gm) and L' (3.78 /zm) 
data have photometric errors smaller than ~0. 1 mag (Bonnefoy 
et al. 2011; Currie et al. 2011b). Other high signal-to-noise 
(S/N) detections such as at M' were obtained without reliable 
flux calibration (Currie et al. 2011b) or with additional, large 
photometric uncertainties due to processing (Bonnefoy et al. 
2013). As a result, the best-fit models admit a wide range of 
temperatures, surface gravities, and cloud structures (e.g., Currie 
et al. 20 1 lb). Thus, new, higher S /N, precise, and flux-calibrated 
photometry at 1-5 /im should provide a clearer picture of 
the clouds, chemistry, temperature, and gravity of ft Pic b. 
Moreover, new near-to-mid IR data may identify distinguishing 
characteristics of /I Pic b’s atmosphere, much like clouds and 
non-equilibrium carbon chemistry for HR 8799 bcde (Currie 
et al. 2011a; Galicher et al. 2011; Skemer et al. 2012; Konopacky 
et al. 2013). 

In this study, we present new 1.5-5 /im observations of 
ft Pic b obtained with NaCo on the VLT and the Near-Infrared 
Coronagraphic Imager (NICI) on Gemini-South. We extract the 
first detection in the 3.09 /zm water-ice filter, the first high S/N, 
well-calibrated H, [4.05], and M' detections, and higher S/N 
detections at K s and L' . We add re-reduced ft Pic data obtained 
in the / (1.25 /zm) and H (1.65 /zm) bands (first presented in 
Bonnefoy et al. 2013) to our new data and we recover ft Pic b 
at a slightly higher S /N and derive its photometry with smaller 
errors. 

We compare the colors derived from broadband photometry 
to those of field substellar objects with a range of spectral types 
to assess whether ft Pic b’s colors appear anomalous/redder than 
the field sequence, as is the case for planets around HR 8799 and 
k And, planet-mass companions like 2M 1207 B, GSC 06214 B, 
and 1RXJ 1609 B (Chauvin et al. 2004; Ireland et al. 2011; 
Lafreniere et al. 2008), and other substellar objects like Luhman 
16B (Luhman et al. 2013). We use atmosphere modeling to 
constrain the range of temperatures, surface gravities, and cloud 
structures plausible for the planet. While previous studies have 
shown the importance of clouds and non-equilibrium carbon 
chemistry in fitting the spectra/photometry of directly imaged 
planets (Bowler et al. 2010; Currie et al. 2011a; Madhusudhan 
et al. 20 1 1 ; Galicher et al. 20 1 1 ; Skemer et al. 2012; Konopacky 
et al. 2013), here the assumed sizes of dust particles entrained 
in the clouds play a critical role. 


2. OBSERVATIONS AND DATA REDUCTION 
2.1. VLT/NaCo Data and Basic Processing 

We observed ft Pic under photometric conditions on 2012 De- 
cember 14 to 17 with the NaCo instrument (Rousset et al. 2003) 
on the VLT UT4/Yepun at Paranal Observatory (Program ID 
090.C-0396). All data were taken in pupil-tracking/angular dif- 
ferential imaging (Marois et al. 2006) and data cube mode. 
Table 1 summarizes the basic properties of these observa- 
tions. Our full complement of data during the run includes 
imaging at 1.04 /zm, 2.12 /tm, K s / 2.18 /zm, 2.32 /zm, 3.74 /zm, 
L'/3.7 8 /zm, Br-a/4.05 /zm, and M' . Here, we focus only on 
the L ' , [4.05], and M' data, deferring the rest to a later study. 
Each observation was centered on ft Pic’s transit for a total field 
rotation of ~50°-70° and total observing times ranging between 
~30 minutes and 59 minutes. 

Along with these new observations, we re-reduce /-band and 
H- band data first presented in Bonnefoy et al. (2013) and taken 
on 2011 December 16 and 2012 January 11, respectively. The 
saturated /-band science images are bracketed by two sequences 
of unsaturated images obtained in a neutral density filter for flux 
calibration. While there were additional frames taken but not 
analyzed in Bonnefoy et al., we found these to be of significantly 
poorer quality and thus do not consider them here. In total, the 
/-band data we consider cover 40 minutes of integration time 
and ~23° of field rotation. The //-hand data cover ~92 minutes 
of integration time and ~36° of field rotation. 

Basic NaCo image processing steps were performed as in 
Currie et al. (2010, 2011b, 2012c) and Rodigas et al. (2012). 
The thermal IR data at L' and [4.05] (M r ) were obtained in 
a dither pattern with offsets every two images (one image) 
to remove the sky background. As all data were obtained in 
data cube mode, we increased our point-spread function (PSF) 
quality by realigning each individual exposure in the cube to 
a common center position and clipping out frames with low 
encircled energy (i.e., those with a core/halo ratio < max(core/ 
halo) — 3 x <r(core-to-halo ratio)). 

2.2. Gemini / NICI Data and Basic Processing 

We obtained Gemini imaging for ft Pic b using the NICI on 
2012 December 23 and 26 in the HiO filter (a„ = 3.09 /zm) and 
on 2013 January 9 in the Hand K s filters (dual-channel imaging), 
both under photometric conditions (Program GS-2012B-Q-40). 
These observations were also executed in angular differential 
imaging (ADI) mode. For the H 2 O data, we dithered each 38 s 
exposure for sky subtraction for a total of ~38 minutes of 
integration time over a field rotation of ~30°. For the H/K s data, 
we placed the star behind the r = 0ft22 partially transmissive 
coronagraphic mask to suppress the stellar halo. Here, we took 
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Table 2 

/! Pic b Detections and Photometry 


UT Date 

Telescope /Instrument 

Filter 

Wavelength 
O m) 

S/N 

Apparent Magnitude 

Absolute Magnitude 

2011 Dec 16 

VLT/NaCo 

J 

1.25 

9.2 

14.11 ±0.21 

12.68 ±0.21 

2012 Jan 1 

VLT/NaCo 

H 

1.65 

30 

13.32 ±0.14 

11.89 ±0.14 

2013 Jan 9 

Gemini/NICI 

H 

1.65 

6.4 

13.25 ±0.18 

11.82 ± 0.18 

2013 Jan 9 

Gemini/NICI 

K s 

2.16 

10 

12.47 ±0.13 

11.04 ±0.13 

2012 Dec 23 

Gemini/NICI 

[3.09] 

3.09 

4.6 

11.71 ±0.27 

10.28 ± 0.27 

2012 Dec 26 

Gemini/NICI 

[3.09] 

3.09 

11 



2012 Dec 16 

VLT/NaCo 

L' 

3.8 

40 

11.24 ±0.08 

9.81 ±0.08 

2012 Dec 16 

VLT/NaCo 

[4.05] 

4.05 

20 

11.04 ±0.08 

9.61 ±0.08 

2012 Dec 15 

VLT/NaCo 

M' 

4.78 

22 

10.96 ±0.13 

9.54 ±0.13 


shorter exposures of fi Pic (f int ~ 11.4 s) to better identify 
and filter out frames with bad AO correction. Our observing 
sequence consists of ~22 minutes of usable data centered on a 
transit with a field rotation of ~41°. 

The basic image processing follows the steps described above 
for the NaCo data. The PSF halo was saturated out to r ~ 
0'.'32-0'.'36 in the H band during most of the observations and 
our sequence suffered periodic seeing bubbles that saturated 
the halo out to angular separations overlapping with the fi Pic b 
PSF. Thus, we focus on reducing only those H- band frames with 
less severe halo saturation (r sat < Of'. 36). The K s observations, 
obtained at a higher Strehl ratio, never suffered halo saturation. 
The first of the two FLO sets suffered from severe periodic seeing 
bubbles and thus generally poor AO performance. We identify 
and remove from analysis frames whose halo flux exceeded 
/’n, m +3fT, where F mm is the minimum flux within an aperture 
covering /I Pic b and a is the dispersion in this flux (about 
10%-25% of the frames, depending on the data set in question). 

2.3. PSF Subtraction 

To remove the noisy stellar halo and reveal /I Pic b, we process 
the data with our “adaptive” locally optimized combination of 
images (A-LOCI) pipeline (Currie et al. 2012a, 2012b; T. Currie 
201 3, in preparation). This approach adopts the LOCI formalism 
(Lafreniere et al. 2007b), where we perform PSF subtraction in 
small annular regions (the “subtraction zone”) at a time over 
each image. We use previously described A-LOCI components, 
including “subtraction zone centering” (Currie et al. 2012b), 
“speckle filtering” to identify and remove images with noise 
structure poorly correlated with that from the science image we 
are wanting to subtract (Currie et al. 2012b), and a moving pixel 
mask to increase point source throughput and normalize it as 
a function of azimuthal angle (Currie et al. 2012a). We do not 
consider a PSF reference library (Currie et al. 20 1 2a) since /I Pic 
is our only target. 

We incorporate a component different from but complemen- 
tary to our “speckle filtering” into A-LOCI, as recently utilized 
in Currie et al. (2012a). We use singular value decomposition 
(SVD) to limit the number of images used in a given annu- 
lar region (i.e., for a given optimization zone) to construct and 
subtract a reference. Briefly, in the (A-)LOCl formalism, a ma- 
trix inversion yields the set of coefficients c k applied to each 
image making up the reference “image”: c k = A 1 1 ). Flere, 
A is the covariance matrix and b is a column matrix defined 
from i pixels in the “optimization zones” of the j th reference 
image section O' and the science image, 0 T : /;, = ]C () J . Of 

i 

(see Lafreniere et al. 2007b). In the previous versions of our 


codes, we use a simple double-precision matrix inversion to in- 
vert the covariance matrix and then solve for c k after multiplying 
by b. 

In this work, we instead use SVD to rewrite A as UEV r 
such that A -1 = VX _ 1 U r , where the T superscript stands for 
the transpose of the matrix. Prior to inversion, we truncate 
the number of singular values at a predefined cutoff, svd\{ m . 
This eigenvalue truncation is very similar to and functions 
the same as the truncation of principal components, /V pca , in 
the Karhunen-Loeve Image Projection (KLIP; Soummer et al. 
2012) and has been successfully incorporated before (Marois 
et al. 2010a). We found that both speckle filtering and SVD 
truncation within our formalism can yield significant contrast 
gains over LOCI and KLIP/Principal Component Analysis 
(PCA), although in this study at the angular separation of fi Pic b 
(«0'.'45), the gains over LOCI are typically about a factor of 1 .5, 
albeit with substantially higher throughput. 1 1 

2.4. Planet Detections 

Figures 1-3 show reduced NaCo and NICI images of fi Pic. 
We detect fi Pic b in all data sets (summarized in Table 2). To 
compute the S/N for fi Pic b, we determine the dispersion, er, 
in pixel values of our final image convolved with a Gaussian 
along a ring with width of 1 full-width at half maximum 
(FWHM) at the same angular separation as fi Pic b but excluding 
the planet (e.g., Thalmann et al. 2009) and we average the 
S/N/pixel over the aperture area. For the Gemini-NICI FI, K s , 
and two [3.1] data sets, the S/Ns are thus 6.4, 11, 4.6, and 
10, respectively. For the J- and H - band NaCo data previously 
presented in Bonnefoy et al. (2013), we achieve S/N ~ 9 and 
S/N ~ 30, respectively. Generally speaking, our 3.8-5 /im 
NaCo data are deeper than the near-IR NaCo data and especially 
the near-IR NICI data, where we detect fi Pic b at S/N = 40 in 
L' and 22 at M', roughly a factor of two higher than previously 
reported (Currie et al. 201 lb; Bonnefoy et al. 2013). These gains 
are due to fi Pic b now being at a wider projected separation 
( L ’ ) or post-processing and slightly better observing conditions 
( M '). The high S /N detections obtained with NaCo also leverage 
on recent engineering upgrades that substantially improved the 
instrument’s image quality and the stability of its PSF (Girard 
et al. 2012). 


11 Recently, Amara & Quanz (2012) claimed a contrast gain of ~5x over 
LOCI using PCA. However, optimal set-ups even within a given formalism 
like LOCI or PCA/KLIP are very data set-specific (cf. Lafreniere et al. 2007b; 
Currie et al. 2012a, 2012b). With LOCI, we obtained roughly equivalent S/Ns 
for ft Pic b during the same observing run but on a night with poorer observing 
conditions (2009 December 29) than a test data set (Currie et al. 201 lb). 
Implementing some A-LOCI filtering and pixel masking yields S/N ~ 30-35. 
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Figure 1 . Processed Gemini-NICI images obtained in the //-band (top left), the /(,-band (top right), and the H 2 O filter centered at ~3.1 // m (bottom panels; 2012 
December 23 data on the left and the 2012 December 26 data on the right). For clarity, we mask the region interior to ~0'.'4 and convolve the image with a Gaussian 
equal to the image FWHM. The planet /J Pic b is in the lower-right at a P.A. of ~210° and a separation of ~0'.'46. The color scale is set such that over the planet’s 
FWHM the pixel color is white. 

(A color version of this figure is available in the online journal.) 



Figure 2. Processed VLT /NaCo ./-(left) and //-band (right) images presented in the same manner as Figure 1 . Owing to very good speckle suppression, the //-band 
data’s effective inner working angle beyond which we are sensitive to fi Pic b-brightness companions is significantly smaller than for the /-band data and the preceding 
NICI images (i'iwa = 072). 

(A color version of this figure is available in the online journal.) 


The optimal A-LOCI algorithm parameters vary significantly 
from data set to dataset. The rotation gap (AP.A., in units of the 
image FWHM) criterion used to produce most of the images is S 
~ 0.6-0.65 , although it is significantly larger for the J and H data 
sets ( 8 — 0.75-0.95). Generally speaking, the optimization areas 
we use Na are significantly smaller (N A = 50-150) than those 
typically adopted (i.e., N A = 300; Lafreniere et al. 2007b). We 
speculate that the pixel masking component of A-LOCI drives 
the optimal N A settings toward these smaller values since the 


planet flux (ostensibly within the subtraction zone) no longer 
significantly biases the coefficient determinations to the point of 
reducing the planet’s S/N. Filtering parameters r C0IT and 
likewise vary wildly from r co „ = 0 and svdum = 2.5 x 10~ 7 at J 
to r corr = 0.9 for the NICI //-band data or .v w/| ml = 2.5 x I 0~ 2 
for the M' NaCo data. 

Although the many algorithm free parameters make finding 
an optimal combination difficult and computationally expensive, 
our final image quality is nevertheless extremely sensitive to 
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Figure 3. Processed VLT/NaCo L' (left), [4.05] (middle), and M' (right) images presented as in Figures 1 and 2. 
(A color version of this figure is available in the online journal.) 


some values, in particular svd\\ m and r COIT . As a test, we explored 
other image processing methods — ADI-based classical PSF 
subtraction and LOCI. While A-LOCI always yields deeper 
contrasts, we easily detect ft Pic b in the mid-IR NaCo data 
using any method and only the poorer of the two [3.1] data sets 
requires A-LOCI to yield a better than 4 a detection (i.e., where 
(Jdet = 1.0857/S/N = 0.27 mag). We will present a detailed 
analysis of image processing methods and algorithm parameters 
in an upcoming study (T. Currie 2013, in preparation). 

Adopting the pixel scales listed in Table 1 , ft Pic b is detected 
at an angular separation of r ~ 0'.'46 in each data set. The 
position angle (P.A.) of ft Pic b is consistent with previous values 
(PA. ^ 210°) and in between values for the main disk and the 
warp, intermediate between the results presented in Currie et al. 
(201 lb) and Lagrange et al. (2012a). While the NICI north P.A. 
on the detector is precisely known and determined from facility 
observations, we have not yet used our astrometric standard 
observations to derive the NaCo P.A. offset, which changes every 
time NaCo is removed from the telescope. To dissuade others 
from using the poorly calibrated NaCo data and the precisely 
calibrated data (Lagrange et al. 2012a) together, we reserve a 
detailed determination of ft Pic b’s astrometry and a study of its 
orbit for a future study. We also detect the ft Pic debris disk in 
each new broadband data set and at [4.05] (Figure 4). We will 
analyze the disk’s properties at a later time as well. 

2.5. Planet Photometry 

To derive ft Pic b photometry, we first measured its brightness 
within an aperture roughly equal to the image FWFIM in each 
case, which was known since we either had AO-corrected 
standard star observations (NICI H, K s , and [3.1]), unsaturated 
images of the primary as seen through the coronagraphic mask 
(NICI K s ), unsaturated neutral density filter observations (NaCo 
7, H, and M'), or unsaturated images of the primary (NaCo 
L' and [4.05]). We assessed and corrected for planet throughput 
losses due to processing by comparing the flux of synthetic 
point sources within this aperture implanted into registered 
images at the same angular separation as ft Pic b before and 
after processing. To derive ft Pic b’s throughput and uncertainty 
in the throughput (<7 a tten), we repeat these measurements at 15 
different P.A.s and adopt the clipped mean of the throughput 
as our throughput and standard deviation of this mean as 
its uncertainty. The planet throughput ranges from 0.38 for 
the 7-band data to 0.82 for the [4.05] data and 0.96 for the 
NICI H - band data, even with aggressive algorithm parameters 



Figure 4. Classical ADI reduction of our L' data showing a clear detection of 
the (3 Pic debris disk. The green dot denotes the position of /? Pic b in the disk. 
(A color version of this figure is available in the online journal.) 


(i.e., <5 ~ 0.6), due to the throughput gains yielded by our pixel 
masking and the SVD cutoff. 

For photometric calibration, we followed several different 
approaches. For the NICI data, we used TYC 7594-1689- 
1 and HD 38921 as photometric standards. We were only 
able to obtain photometric calibrations for the first of the 
two [3.1] data sets. For all other data, we used the primary 
star, ft Pic, for flux calibration, adopting the measurements 
listed in Bonnefoy et al. (2013). For the 7 and H NaCo 
data, we used images of the primary as viewed through a 
neutral density filter. For the M' NaCo data, we obtained 
neutral density filter observations and very short exposures. 
While the latter were close to saturation and were probably 
in the non-linear regime, the implied photometry of ft Pic was 
consistent to within the errors. The primary was unsaturated 
in the [4.05] data. Finally, for the L' data, we took 8.372 ms 
unsaturated images of ft Pic for flux calibration. In all cases, 
we again adopt the standard deviation of the clipped mean 
of individual measurements as our photometric calibration 
uncertainty, erg uxca i. To compute the photometric uncertainty 
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Table 3 

Sample Photometric NaCo and NICI Error Budgets 


Telescope /Instrument 

Filter 

Apparent Magnitude 

O/iet 

fatten 

O’fluxcal 

Gemini/NICI 

K s 

12.47 ± 0.13 

0.11 

0.06 

0.03 

Gemini/NICI 

[3.09] 

11.71 ± 0.27 

0.24 

0.11 

0.08 

VLT/NaCo 

L' 

11.24 ±0.08 

0.03 

0.07 

0.03 

VLT/NaCo 

M' 

10.96 ±0.13 

0.05 

0.11 

0.06 


Note. The photometric uncertainty from the intrinsic S/N (er^et) scales as 
1.0857/S/N. 


for each data set, we considered the S/N of our detection, 
the uncertainty in the planet throughput, and the uncertainty 
in absolute flux calibration: cr = V cr 2 et + tr 2 ten + ff fl 2 uxcal . 

Table 2 reports our photometry and Table 3 lists sample 
error budgets for the two NICI photometric measurements and 
the two NaCo measurements. The relative contributions from 
each source of photometric uncertainty to the total uncertainty 
are representative of our combined data set. For the [3.09] 
data, residual speckle noise/sky fluctuations greatly limit the 
planet’s S /N and thus crdet is the primary source of photometric 
uncertainty. For the K s data, the intrinsic S /N and the two other 
sources of photometric uncertainty contribute in more equal 
proportions. The L' and M' data error budgets are characteristic 
of most of our other data, where the photometric uncertainty 
is primarily due to the absolute photometric calibration and 
throughput. With the exception of the [3.09] NICI data, the 
intrinsic S/N of the detection does not dominate the error 
budget. For the best quality (mid-IR NaCo) data, the throughput 
uncertainty was small (~5%') and was never any larger than 
15% (/band data) in any data set. 12 

In general, we find fair agreement with previously published 
photometry, where our measurements are usually consistent 
within the photometric errors with those reported previously 
(e.g., m H = 13.32 ± 0.14 and 13.25 ± 0.18 versus 13.5 ± 0.2 
in Bonnefoy et al. 2013). Our L' photometry is more consistent 
with Currie et al.’s measurement of mu— 9.73 ± 0.06 than with 
that listed in Bonnefoy et al. (2013) (mu— 9.5 ± 0.2), although it 
is nearly identical to that derived for some ft Pic b data sets listed 
in Lagrange et al. (2010). Our [4.05] photometry implies that 
f J > Pic b is ~15%-20% brighter there than previously assumed 
(Quanz et al. 2010) and may have a slightly red L'— [4.05] 
color. The major difference from previous studies, however, is 
that our photometric errors are consistently much smaller. For 
example, the uncertainty on the [4.05] photometry is reduced to 
0.08 mag from 0.23 mag due both to higher S/N detections and 
lower uncertainty in our derived photometry (i.e., throughput 
corrections). The NICI photometry is also substantially less 
uncertain than given in Boccaletti et al. (2013) because f Pic b is 
not occulted by the focal plane mask. These lower uncertainties 


12 In principle, tuning the algorithm parameters to maximize the S/N of 
f! Pic b could introduce additional photometric uncertainties if the planet is in 
significant residual speckle contamination. In such a case, the algorithm 
parameters maximizing the S/N could instead be the set that maximizes the 
residual speckle contamination within the planet aperture while minimizing it 
elsewhere, especially since the pixel masking technique normalizes the point 
source throughput but not the noise as a function of azimuthal angle. However, 
we do not find substantial differences in the derived photometry if we adopt a 
default set of algorithm parameters. Furthermore, the parameters maximizing 
the S/N are never the ones maximizing the planet throughput and our tuning is 
not just finding the parameter set making pixels within the planet aperture the 
“noisiest.” Adopting slightly different parameters from the “optimized” case 
yields nearly identical photometry. Moreover, residual speckle contamination 
in most data sets is extremely low, the intrinsic S/N is limited by sky 
background fluctuations in addition to speckles for the mid-IR data. 


should allow more robust comparisons between ft Pic b and 
other substellar objects and, from modeling, more precise limits 
on the best-fitting planet atmosphere properties. 

3. EMPIRICAL COMPARISONS WITH f Pic b 

Our new data allow us to compare the spectral energy dis- 
tribution (SED) of ft Pic b to the SEDs of many field L/T- 
type brown dwarfs, as well as directly imaged, low-surface 
gravity, low-mass brown dwarf companions and directly im- 
aged planets. Our goal here is to place f Pic b within the 
general L/T-type spectral sequence, identify departures from 
this sequence such as those seen for low surface gravity ob- 
jects like HR 8799 bcde, and identify the substellar object(s) 
with the best-matched SED. Some bona fide directly imaged 
planets like HR 8799 bcde and at least some of the lowest- 
mass brown dwarfs like 2M 1207 B appear redder/cloudier 
than their field dwarf counterparts at comparable temperatures 
(7/ re ~ 900-1100 K). However, it is unclear whether hotter 
imaged exoplanets appear different from their (already cloudy) 
field L dwarf counterparts and f Pic b provides a test of any such 
differences. We will use our comparisons with the L/T dwarf 
sequence and the SEDs of other substellar objects to inform our 
atmosphere model comparisons later to derive planet physical 
parameters (e.g., 7/ff and log(g)). 

3.1. Infrared Colors of fi Pic b 

To compare the near-to-mid IR properties of ft Pic b with 
those of other cool, substellar objects, we primarily use the 
sample of L/T dwarfs compiled by Leggett et al. (2010), which 
includes field dwarf spectral classes between ~M7 and T5, 
corresponding to a range of temperatures between ~2500 K and 
700 K. To explore how the ft Pic b SED compares to those with 
other directly imaged planets/planet candidates and very low- 
mass brown dwarf companions within this temperature range, 
we include objects listed in Table 4. These include the directly 
imaged planets around HR 8799 (Marois et al. 2008, 2010b; 
Currie et al. 2011a) and the directly imaged planet candidate 
around k And (Carson et al. 2013). Additionally, we include 
high mass ratio brown dwarf companions with masses less 
than the deuterium-burning limit (~ 13—14 Mj) and higher- 
mass companions whose youth likely favors a lower surface 
gravity than for field brown dwarfs, a difference that affects the 
objects’ spectra (e.g., Luhman et al. 2007). Among these objects 
are 1RXJ 1609B, AB Pic B, and Luhman 16 B (Lafreniere 
et al. 2008; Chauvin et al. 2005; Luhman et al. 2013). Table 5 
compiles photometry for all of these low surface gravity objects. 

Figure 5 compares the IR colors of f Pic b (dark blue 
diamonds) with those of field M dwarfs (small black dots), field 
L0-L5 dwarfs (gray dots), field L5.1-L9 dwarfs (asterisks), T 
dwarfs (small, light gray dots), and planets/low-mass young 
brown dwarfs (light blue squares). The J — H / H — K s colors 
of f Pic b appear slightly blue in J — H and red in H — K s 
compared with field L0-L5 dwarfs, although the difference 
here is not as large as was found by Bonnefoy et al. (2013). 
Other young substellar objects appear to have similar near- 
IR colors, in particular k And b, GSC 06214 B, USco CTIO 
108B, 2M 1207 A, and Luhman 16 B, whose spectral types 
range between M8 and TO. 5. 

The mid-IR colors of f Pic b (top-right and bottom panels) 
show a more complicated situation. In / — K s / K s — L' and 
H — K s /K s — L ' , ft Pic b lies along the field L/T dwarf locus 
with colors in between those for L0-L5 and L5.1-L9 dwarfs, 
overlapping in color with k And b, 1RXJ 1609B, GSC 06214B, 
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Table 4 

Young Directly Imaged Planets and Very Low Mass Brown Dwarfs Used for Comparison 


Companion 

D 

(pc) 

Assoc. 

Age 

(Myr) 

ST(Primary) 

ST(Companion) 

Sep. 

(Au) 

Mass 

(Mj) 

References 

Planets and Planet Candidates 








HR 8799 b 

39.4 ± 1 

Columba 

30 

A5 

?? 

67.5-70.8 

4-5 

1,2, 3, 4, 5 

HR 8799 c 

39.4 ± 1 

Columba 

30 

A5 

L/T? 

42.1^14.4 

-7 

1,2,3, 4, 5 

HR 8799 d 

39.4 ± 1 

Columba 

30 

A5 

L/T? 

26.4-28.1 

-7 

1,2,3, 4, 5 

HR 8799 e 

39.4 ± 1 

Columba 

30 

A5 

L/T? 

-15 

-7 

2,4,6 

k And b 

51.6 ±0.5 

Columba 

30 

B9TV 

L2-L8? 

55 ±2 

11.8-14.8 

7 

Low-mass Brown Dwarfs 








1RXJ 1609 B 

145 ± 14 

us 

5-10 

K7 

L4 ± 2 

330 

6-12 

8,9 

GSC 06214 B 

145 ± 14 

US 

5-10 

K7 

L0± 1 

320 ± 30 

14 ±2 

10,11 

USco CTIO 108 B 

145 ± 14 

US 

5-10 

M7 

M9.5 

670 ± 64 

6-16 

12 

HIP 78530 B 

156.7 ± 13.0 

US 

5-10 

B9V 

M8 ± 1 

710 ±60 

22 ±4 

11,13 

2M 1207 B 

52.4 ± 1.1 

TWA 

—8 

M8 

?? 

40.8 ± 9 

8 ± 2 

14 

2M 1207 A 

52.5 ± 1.1 

TWA 

—8 

M8 

M8 

40.8 ± 9 

24 ±6 

14 

TWA5B 

44.4 ± 4 

TWA 

—8 

M2Ve 

M8-M8.5 

-98 

-20 

15 

HR 7329 B 

47.7 ± 1.5 

ft Pic 

12 

A0 

M7.5 

200 ±7 

26 ±4 

16 

PZ Tel B 

51.5 ± 2.6 

ft Pic 

12 

K0 

M7 ± 2 

17.9 ± 0.9 

36 ±6 

17 

2M0103AB B 

47.2 ±3.1 

Tuc-Hor? 

30 

M5/M6 

L? 

84 

12-14 

18 

AB Pic B 

45.5 ± 1.8 

Carina 

30 

KIVe 

L0± 1 

248 ± 10 

13-14 

19, 20 

Luhman 16 B 

2.02 ±0.15 

Argus? 

40? 

L7.5 

TO. 5 

3.12 ±0.25 

40-65 

21,22 

Luhman 16 A 

2.02 ±0.15 

Argus? 

40? 

L7.5 

L7.5 

3.12 ±0.25 

40-65 

21,22 

CD-35 2722 B 

21.3 ± 1.4 

AB Dor 

-100 

MlVe 

L4± 1 

67.4 ± 4 

31 ±8 

23 


References. (1) Marois et al. 2008; (2) Currie et al. 2011a; (3) Galicher et al. 2011; (4) Skemer et al. 2012; (5) Currie et al. 2012b; (6) Marois et al. 2010b; (7) 
Carson et al. 2013; (8) Lafreniere et al. 2008; (9) Lafreniere et al. 2010; (10) Ireland et al. 2011; (11) Bailey et al. 2013; (12) Bejar et al. 2008; (13) Lafreniere 
et al. 2010; (14) Chauvin et al. 2004; (15) Lowrance et al. 1999; (16) Lowrance et al. 2000; (17) Biller et al. 2010; (18) Delorme et al. 2013; (19) Chauvin 
et al. 2005; (20) Bonnefoy et al. 2010; (21) Luhman et al. 2013; (22) Burgasser et al. 2013; (23) Wahhaj et al. 2011. 


Table 5 

Photometry for Young Directly Imaged Planets and Very Low Mass Brown Dwarfs 


Companion 

J 

H 

K s 

[3.09] 

L' 

M' 

Z^Pic b 

C.L. 

Planets and Planet Candidates 








HR 8799 b 

16.52 ±0.14 

15.08 ±0.17 

14.05 ± 0.08 


12.68 ±0.12 

13.07 ± 0.30 

52.8 

-1 

HR 8799 c 

14.65 ±0.17 

14.18 ± 0.17 

13.13 ±0.08 


11.83 ± 0.07 

12.05 ±0.14 

6.098 

0.893 

HR 8799 d 

15.26 ±0.43 

14.23 ± 0.22 

13.11 ±0.08 


11.50 ±0.12 

11.67 ±0.35 

8.351 

0.961 

HR 8799 e 


13.88 ± 0.20 

12.89 ± 0.26 


11.61 ± 0.12 

>10.09 



k And b 

12.7 ± 0.30 

11.7 ±0.20 

11.0 ±0.4 


9.54 ± 0.09 


0.946 

0.186 

Low-mass Brown Dwarfs 








1RXJ 1609 B 

12.09 ±0.12 

11.06 ±0.07 

10.38 ± 0.05 

9.84 ±0.21 

8.99 ± 0.30 


1.369 

0.287 

GSC 06214 B 

10.43 ± 0.04 

9.74 ± 0.04 

9.14 ±0.05 

8.60 ± 0.08 

7.94 ± 0.07 

7.94 ± 0.30 

7.001 

0.864 

USco CTIO 108 B 

10.72 ± 0.09 

9.94 ± 0.08 

9.30 ±0.11 






HIP 78530 B 

9.25 ± 0.05 

8.58 ± 0.04 

8.36 ± 0.04 


7.99 ± 0.06 


88.086 

-1 

2M 1207 B 

16.40 ± 0.2 

14.49 ±0.21 

13.31 ±0.11 


11.68 ±0.14 


20.601 

-1 

2M 1207 A 

9.35 ± 0.03 

8.74 ± 0.03 

8.30 ± 0.03 


7.73 ± 0.10 


48.044 

-1 

TWA 5 B 

9.1 ±0.2 

8.65 ± 0.06 

7.91 ±0.2 






HR 7329 B 

8.64 ±0.19 

8.33 ± 0.1 

8.18 ±0.1 


7.69 ±0.1 


59.455 

-1 

PZ Tel B 

8.70 ±0.18 

8.31 ± 0.15 

7.86 ±0.19 






2M0103AB B 

12.1 ±0.3 

10.9 ± 0.2 

10.3 ± 0.2 


9.3 ± 0.1 


2.666 

0.736 

AB Pic B 

12.80 ±0.10 

11.31 ± 0.10 

10.76 ± 0.08 


9.9 ±0.1 


11.231 

0.996 

Luhman 16 B 

14.69 ± 0.04 

13.86 ± 0.04 

13.20 ± 0.09 






Luhman 16 A 

15.00 ±0.04 

13.84 ± 0.04 

12.91 ±0.04 






CD-35 2722 B 

11.99 ±0.18 

1 1.14 ± 0.19 

10.37 ±0.16 







Note. We only quantitatively compare the photometry between ft Pic b and those objects with full JHK S L' photometry. 


HR 8799 d, and 2M 1207 B. Compared with the few field 
L/T dwarfs from the Leggett et al. sample with M' photometry, 
ft Pic b appears rather red, most similar in K s — M' color to 
GSC 06214 B. 

The color-magnitude diagram positions of ft Pic b (Figure 6) 
better clarify how its near-to-mid SED compares to the field 


L/T dwarf sequence and to very low mass (and gravity?) 
young substellar objects. In general, compared to the field L 
dwarf sequence, ft Pic b appears progressively redder at mid-IR 
wavelengths. Similar to the case for GSC 06214 B (Bailey et al. 
2013), ft Pic b appears overluminous compared to the entire 
L/T dwarf sequence in the mid-IR. 
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Figure 5. Near-to-mid infrared color-color diagrams comparing ft Pic b's colors (blue diamonds) with those of M dwarfs (small dark circles), early L-type brown dwarfs 
(gray circles), late L dwarfs (asterisks), and T dwarfs (light gray dots), from Leggett et al. (2010). We also overplot the positions of young, substellar objects/other 
directly imaged planets (turquoise squares). 

(A color version of this figure is available in the online journal.) 


3.2. Comparisons with SEDs of Other Substellar Objects 

To further explore how the SED of j5 Pic b agrees with/departs 
from the field L/T dwarf sequence and other young substellar 
objects, we first compare its photometry to spectra from the 
SpeX library (Cushing et al. 2005; Rayner et al. 2009) of brown 
dwarfs with data overlapping with our narrowband mid-IR filters 
([3.09] and [4.05]) spanning spectral classes between LI and L5: 
2MASS J14392836+1929149 (LI), Kelu-lAB (L2), 2MASS 
J 15065441+1 321060 (L3), and 2MASS J15074769- 1627386 
(L5). To compare the f Pic b photometry with cooler L 
dwarfs, we add combined IRTF/SpeX and Subaru/IRCS spectra 
from 1 to 4.1 /xm for 2MASS J0825 1968+21 15521 (L7.5) and 
DENIS-P J025503. 3— 470049 (L8) (Cushing et al. 2008). 
Finally, we add spectra for the low surface gravity L4.5 dwarf, 
2MASS J2224438 1—0158521 (Cushing et al. 2008). To high- 
light differences between f> Pic b and these L dwarfs, we scale 
the flux densities for each of these standards to match f> Pic b at 
~2.15 p , m ( K s band). 

To convert our photometry derived in magnitudes to flux 
density units, we use the zeropoint fluxes listed in Table 6. 
The JHK S and i'M'(4.78 /xm) zeropoints are from Cohen et al. 
(2003) and Tokunaga & Vacca (2005), respectively. We base the 
other zeropoints off of Rieke et al. (2008), although alternate 
sources (e.g., Cohen et al. 1995) yield nearly identical values. 


Table 6 

Adopted Flux Density Zeropoints 


Filter 

(/xm) 

F v ,o 

(Jy) 

J 

1.25 

1594 

H 

1.65 

1024 

K s 

2.15 

666.20 

[3.09] 

3.09 

356 

L' 

3.78 

248 

[4.05] 

4.05 

207 

M’ 

4.78 

154 


Because the overlap in wavelengths between ft Pic and these 
objects is not uniform, we do not perform a rigorous fit between 
the two, finding the scaling factor that minimizes the x 2 value 
defined from the planet flux density, the comparison object 
flux density, and the photometric errors in both. Rather, we 
focus on a simple first-order comparison between f> Pic b and 
the comparison objects to motivate the detailed atmospheric 
modeling later in Section 4. 

Figure 7 (left panel) compares the photometry of f> Pic b with 
spectra for field L1-L5 dwarfs. While the LI standard slightly 
overpredicts the flux density in the J band, the other three 
early/mid L standards match the f Pic b near-lR SED quite 
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Figure 6. Color-magnitude diagrams comparing ft Pic b with field M, L, and T dwarfs and young, substellar objects. Symbols are the same as in Figure 5. 

(A color version of this figure is available in the online journal.) 



Figure 7. Photometric data for ft Pic b compared with L dwarf standard spectra between L2 and L8, as well as a low surface gravity L4 dwarf from Cushing et al. 
(2005, 2008). We scale ft Pic b’s flux density in the A" 5 -band to the band-integrated flux density of the standards at 2.15 /xm. With the possible exception of the low 
surface gravity L4 dwarf, none of these standards provide a good match to the measurements of ft Pic b. We identify the passbands along the bottom of the plot 
(left-hand panel). 

(A color version of this figure is available in the online journal.) 


well, indicating a “near-IR spectral type” of ~L2-L5. The L7.5 
and L8 standards also produce reasonable matches, although 
they tend to underpredict the brightness in the J band (right 
panel). 


However, all standards have difficulty matching the ft Pic b 
SED from 3-4 /rm. In particular, the ft Pic b flux density from 
~3 to ~5 ji m is nearly constant, whereas it rises through 4 // m 
and then steeply drops in all six standards depicted here. When 
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Figure 8. Comparisons between the ft Pic b SED and the closest-matching 
substellar objects with JHK S L' photometry. Quantitatively, 1RXJ 1609 B and 
k And b provide the closest matches, although it is, as yet, unclear whether 
any known substellar object fully reproduces ft Pic b’s SED at all measured 
wavelengths. 

(A color version of this figure is available in the online journal.) 


one focuses on only the p Pic b photometry over 3.8-4. 1 /im, 
the “mid-IR spectral type” is hard to define; the low surface 
gravity L4.5 dwarf bears the greatest resemblance, although we 
fail to identify good matches at all wavelengths with any of 
our spectral templates, where the 3.1 /xm, L', and [4.05] data 
points are the most problematic. While none of our standards 
have measurements fully overlapping with the M' filter, the 
flux densities at 5.1 /xm indicate that the standards may have a 
very hard time simultaneously reproducing our measurements 
in all four filters between 3 and 5 /xm. Although non-equilibrium 
carbon chemistry can flatten the spectra of low surface gravity 
L/T dwarfs (Skemer et al. 2012), its effect is to weaken the 
methane absorption trough at ~3.3 /xm and suppress emission 
at ~5 /xm. Thus, it is unclear whether this effect can explain 
the enhanced emission at ~3.1 /xm (mostly outside of the CH4 
absorption feature to begin with) and 5 /xm. 

To understand whether ( J > Pic b’s SED is unique even amongst 
other very low-mass substellar objects, we compare our pho- 
tometry with that of the companions listed in Table 4 that have 
photometry from 1 /xm through ~4-5 /xm: HR 8799 bed, k 
And b, 1RXJ 1609 B, GSC 06214B, HIP 78530 B, 2M 1207A/ 

B, HR 7329B, and AB Pic. Two objects— 1RXJ 1609 B and 
GSC 06214B — have 3.1 /xm photometry, k And b has [4.05] 
from data obtained by T. C. (M[4 05] = 9.45 ± 0.20) (Bonnefoy 
et al. 2013). 

The two far-right columns of Table 5 list the reduced / 2 
and the goodness-of-fit statistics between p Pic b’s JHK S L ’ 
([3.1], [4.05]) photometry, while Figure 8 displays these com- 
parisons for k And b, 1RXJ 1609B, and GSC 06214B, which 
are all thought to be low surface gravity companions with 
T eS ~ 1700 K, 1800 K, and 2200 K (Carson et al. 2013; 
Lafreniere et al. 2010; Bowler et al. 201 1; Bailey et al. 2013), re- 
spectively. Overall, k And b provides the best match to ft Pic b’s 
photometry, requires negligible flux scaling, and is essentially 
the same within the 68% confidence limit (C.L.) (/ 2 = 0.946, 

C. L. = 0.186), although the large photometric uncertainties in 
the near-IR limit the robustness of these conclusions. The com- 
panion to 1RXJ 1609 likewise produces a very good match 
(X 2 = 1.369, C.L. = 0.287), while the slightly more lumi- 


nous (and massive) GSC 062 14B appears to be much bluer and 
(relatively) overluminous in V and M' (or, conversely, over- 
luminous at JHK S ) by ~30%. In comparison, the cooler (7’ c n 
Rs 900-1100 K) exoplanets HR 8799 bed provide far poorer 
matches (x 2 ~ 6-52). 

Still, it is unclear whether any object matches p Pic b’s pho- 
tometry at all wavelengths; both of the objects for which we have 
[3.1] data, GSC 06214B and 1RXJ 1609B, are still slightly un- 
derluminous here. Moreover, the best-matching companions — 
a: And b and 1RXJ 1609B — are still not identical, as the scal- 
ing factors between ft Pic b’s spectrum and these companions’ 
spectra that minimize x 2 are ~0.83 and 0.53, respectively. While 
companions with identical temperatures but radii 10% and 30% 
larger than p Pic b would achieve this scaling, k And b and 
1RXJ 1609B are older and younger, respectively, than p Pic b, 
whereas for a given initial entropy of formation planet radii are 
expected to decrease with time (Spiegel & Burrows 2012). 

In summary, young (low surface gravity?), low-mass objects 
may provide a better match to p Pic b’s photometry than do field 
dwarfs, especially those with temperatures well above 1000 K 
but slightly below 2000 K (k And b, 1RXJ 1609 B). However, 
we fail to find a match (within the error bars) between the 
planet’s photometry spanning the full range of wavelengths for 
which we have data, especially at ~3 /xm. As planet spectra 
depend critically on temperature, surface gravity, clouds, and 
(as we shall see) dust particle sizes, our comparisons imply that 
P Pic b may differ from most young substellar objects in one of 
these respects. Next, we turn to detailed atmospheric modeling 
to identify the set of atmospheric parameters that best fit the 
P Pic b data. 

4. PLANET ATMOSPHERE MODELING 

To further explore the physical properties of p Pic b, we com- 
pare its photometry with planet atmosphere models adopting 
a range of surface gravities, effective temperatures, and cloud 
prescriptions/dust. For a given surface gravity and effective 
temperature, a planet’s emitted spectrum depends primarily on 
the atmosphere’s composition, the structure of its clouds, and 
the sizes of the dust particles of which the clouds are comprised 
(Burrows et al. 2006). For simplicity, we assume solar abun- 
dances except where noted and leave consideration of anoma- 
lous abundances for future work. 

Based on p Pic b’s expected luminosity (log(L p /L Q ) ~ —3.7 
to —4, Lagrange et al. 2010; Bonnefoy et al. 2013) and age, it 
is likely too hot (7 [4 ~ 1400-1800 K) for non-equilibrium 
carbon chemistry to play a dominant role (Hubeny & Burrows 
2007; Galicher et al. 2011). Therefore, our atmosphere models 
primarily differ in their treatment of clouds and the dust particles 
entrained in clouds. For each model, we explore a range of 
surface gravities and effective temperatures. 

4.1. Limiting Cases: The Burrows et al. (2006) E60 and A60 
Models and the AMES-DUSTY Models 

4.1.1. Model Descriptions 

We begin by applying an illustrative collection of previously 
developed atmosphere models to p Pic b. These models will 
produce limiting cases for the planet’s cloud structure and 
typical dust grain size, which we refine in Section 4.2. To 
probe the impact of cloud thickness, we first adopt a (large) 
modal particle size of 60 /xm and consider three different cloud 
models: the standard chemical equilibrium atmosphere thin- 
cloud models from Burrows et al. (2006), which successfully 
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Table 7 

f5 Pic b Atmosphere Modeling Grid 


Model 

Teff 

Range 

log(g) 

R P ( Rj ) 

References 

Limiting Cases 

E60 

1000-1800 

4^1.5 

0.9-2 

i 

AE60 

1000-1700 

4^1.5 

0.9-2 

2 

A60 

1000-1700 

4^1.5 

0.9-2 

1,3 

AMES-DUSTY 

1000-2000 

3. 5^4.5 

0.9-2 

4 

New Models 

A4 

1400-1900 

3.6^1.25 

0.9-2 a 

5 


Note. 

a In our modeling, we perform two runs for the A4 models: one where we fix 
the radius to values adopted in the Burrows et al. (1997) hot-start models and 
one where we allow the radius to freely vary between 0.9 Rj and 2 Rj. 
References. (1) Burrows et al. 2006; (2) Madhusudhan et al. 2011; (3) Currie 
et al. 2011a; (4) Allard et al. 2001; (5) this work. 

reproduce the spectra of field L dwarfs, moderately thick cloud 
models from Madhusudhan et al. (2011), and thick cloud models 
used in Currie et al. (2011a). To investigate the impact of 
particle size, we then apply the AMES-DUSTY models. The 
DUSTY models lack any dust grain sedimentation, such that 
the dust grains are everywhere in the atmosphere, similar to the 
distribution of dust grains entrained in thick clouds. However, 
they adopt far smaller dust grains than do the thick cloud 
models from Madhusudhan et al. (2011) and Currie et al. 
(2011a), where the grains are submicron in size and follow 
the interstellar grain size distribution (Allard et al. 2001). All 
models described here and elsewhere in the paper assume that 
the planet is in hydrostatic and radiative equilibrium. None of the 
models consider irradiation from the star, as irradiation is likely 
unimportant at f Pic b’s orbital separation. Table 7 summarizes 
the range of atmospheric properties that we consider for each 
model. 

The Burrows et al. (2006) E60 thin cloud, large dust particle 
models. As described in Burrows et al. (2006) and later works 
(e.g., Currie et al. 201 la; Madhusudhan et al. 2011), the Model 
E60 case assumes that the clouds are confined to a thin layer 
where the thickness of the flat part of the cloud encompasses 
the condensation points of different species with different 
temperature-pressure point intercepts. Above and below this 
flat portion, the cloud shape function decays as the —6th and 
— 10th power, respectively, so that the clouds have scale heights 
of ~l/7th and 1/1 1th that of the gas. We adopt a modal 
particle size of 60 fim and a particle size distribution drawn 
from terrestrial water clouds (Deirmendjian 1964). We consider 
surface gravities with log(g) = 4 and 4.5 and temperatures with 
a range of T c ff = 1400-1800 K in increments of 100 K. 

The Madhusudhan etal. (2011) AE60 moderately thick cloud, 
large dust particle models. The Model AE60 case, described in 
Madhusudhan et al. (2011), assumes a shallower cloud shape 
function of s u = 1, such that the cloud scale height is half that 
of the gas as a whole. We again adopt a modal particle size 
of 60 /im and the same particle size distribution. We consider 
surface gravities with log(g) = 4 and 4.5 and temperatures 
between 7’ c ff = 1000-1700 K in increments of 100 K. 

The Burrows et al. ( 2006 ) A60 thick cloud, large dust particle 
models. As described in Currie et al. (2011a), the Model A60 
case differs in that it assumes that the clouds extend with a scale 
height that tracks that of the gas as a whole. Below the flat part 
of the cloud, the shape function decays as the —10th power as 


in the E60 and AE60 models (Burrows et al. 2006), although 
deviations from this scaling do not affect the emergent spectrum. 
Here, we consider surface gravities with log(g) = 4 and 4.5 and 
temperatures with a range of T e g — 1000-1700 K in increments 
of 100 K. 

AMES-DUSTY thick cloud, small dust particle limit. The 
AMES-DUSTY atmosphere models (Allard et al. 2001) lever- 
age on the PHOENIX radiative transfer code (Hauschildt & 
Baron 1999) and explore the limiting case where dust grains do 
not sediment/rain out in the atmosphere. Unlike the Burrows 
et al. (2006) models and those considered in later works (e.g., 
Spiegel & Burrows 2012), the AMES-DUSTY models adopt 
an interstellar grain size distribution favoring far smaller dust 
grains with higher opacities. The grains’ higher opacities reduce 
the planet’s radiation at shorter wavelengths. Thus, these mod- 
els have dramatically different near-IR planet spectra from the 
E/A/AE60 type models with larger modal grain sizes even at 
the same temperatures and gravities (cf. Burrows et al. 2006; 
Currie et al. 2011a). Here, we consider AMES-DUSTY mod- 
els with log(g) = 3.5, 4, and 4.5 and T eif = 1000-2000 K 
(AT eff =100 K). 

4.1.2. Fitting Method 

To transform the DUSTY spectra into predicted flux den- 
sity measurements (at 10 pc), we convolve the spectra over 
the filter response functions and scale by a dilution factor 
f — (Rpi an et/10 pc) 2 . We consider a range of planet radii be- 
tween 0.9 Rj and 2 Rj. Likewise, we convolve the E60 and A60 
model spectra over the filter response functions. The E60 mod- 
els adopt a mapping between planet radius and surface gravity/ 
temperature set by the Burrows et al. (1997) planet evolution 
models, as do all other Burrows et al. (2006) and Madhusudhan 
et al. (2011) models. To explore departures from these models, 
we allow the radius to vary by an additional scale factor ranging 
from 0.7 to 1.7. For most of our grid, this range translates into 
a radius range of 0.9 to 2 Rj. 

Our atmosphere model fitting follows methods in Currie et al. 
(2011a, 2011b), where we quantify the model fits with the y 1 
statistic, 

n 

X~ “ ^ ' (/data.i ^Znodel,;)~/°data,i ■ (I) 

i=0 

We weight each data point equally. Because our photometric 
calibration fully considers uncertainties due to the S/N, the 
processing-induced attenuation, and the absolute photometric 
calibration, we do not set a 0.1 mag floor to a for each data 
point as we have done previously. 

We determine which models are formally consistent with the 
data by comparing the resulting / 2 value to that identifying the 
68% confidence limit and identify those that can clearly be ruled 
out by computed the 95% confidence limit. Note here that these 
limits are significantly more stringent compared with the ones 
we adopted in Currie et al. (2011a). Treating the planet radius 
as a free parameter, we have five degrees of freedom for seven 
data points, leading to x£g% = 5.87 and /| 5% = 1 1.06. 

4.1.3. Results 

Table 8 summarizes our fitting results using the E60, AE60, 
A60, and DUSTY models. Figure 9 displays some of these fitting 
results, where the left-hand panels show the x 2 distributions with 
the 68% and 95% confidence limits indicated by the horizontal 
dashed and dotted lines, respectively. The right-hand panels 
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Table 8 

Model Fitting Results 


Currie et al. 


Model 

Xmin 

log(g), T tff , R(Rj) (for xiL) 

log(g), r eff , R(Rj) (68%) 

log(g), Tefl, R(Rj) (95%) 

E60 

53.80 

4.0, 1400 K, 1.79 



AE60 

37.70 

4.5, 1400 K, 1.96 



A60 

19.99 

4.5, 1400 K, 2.05 



AMES-DUSTY 

7.14 

3.5, 1700 K, 1.35 


3.5-4, 1700 K, 1.32-1.35 

A4 (fixed radius) 

5.85 

3.8, 1600 K, 1.65 

3.8, 1600 K, 1.65 

3.6, 1500-1550 K, 1.79-1.80 
3.8, 1550-1625 K, 1.65 
4.0, 1650-1700 K, 1.54 

A4 (scale) 

5.82 

4-4.25, 1600 K, 1.65 

3.6^1.25, 1600 K, 1.64-1.66 

3.6-4.25, 1500-1750 K, 1.44-1.82 

A4 (1% CH 4 , scale) 

5.13 

4, 1575 K, 1.71 

4, 1575-1650 K, 1.59-1.71 

4, 1575-1650 K, 1.59-1.71 


Notes. The x 2 values quoted here do not refer to / 2 per degree of freedom. The columns for A4 (fixed radius) do not have a range in radius in some 
columns because only one model (with a fixed radius) identifies the x 1 minimum. 


and middle-left panel show the best-fitting models for each 
atmosphere prescription. A successful model must match three 
key properties of the observed SED: (1) At 3-5 /im, the SED 
is relatively flat; (2) at 1-3 /tm, the spectral slope is relatively 
shallow; and (3) the overall normalization of the 3-5 /im flux 
relative to the 1-3 /im flux must match the data. 

For the E60, AE60, and A60 models, we find y 2 minima 
at log(g) = 4-4.5 and 7/® = 1400 K in each case with 
radius scaling factors, the constant that we multiply the nominal 
Burrows et al. planet radii by, between 1.185 and 1.680. For the 
Burrows et al. (1997) evolutionary models, these scaling factors 
imply planet radii between ~1.8 and 2 Rj, at the upper extrema 
of our grid in radius. 

Figure 9 illustrates the impact of changing cloud models on 
the SED, given a fixed grain size. The best-fit temperature does 
not vary dramatically because, roughly speaking, the relative 
fluxes at 1-3 /im and 3-5 /im are determined by the SED’s 
blackbody envelope. However, cloud thickness dramatically 
affects the depths of absorption bands superimposed on that 
envelope. Atmosphere models presented here do not feature 
temperature inversions. As such, high opacity molecular lines 
have low flux densities because they originate at high altitudes 
where the temperature is low. When clouds are thin, an optical 
depth of unity is achieved at very different altitudes in and 
outside of absorption bands such as those at 3.3 jim (methane) 
and 4.5 /im (primarily CO) and the bands appear deep. 

For a fixed observed effective temperature, thicker clouds 
translate into hotter temperature profiles (i.e., at a given pressure 
in the atmosphere, the temperature is higher) (e.g., Madhusud- 
han et al. 201 1). The total Rosseland mean optical depth of the 
atmosphere at a given pressure is higher (Madhusudhan et al. 
2011). As the clouds become thicker, the r = 1 surface also 
is more uniform, such that molecular features wash out and the 
spectrum overall appears flatter and more like a blackbody (Bur- 
rows et al. 2006). Hence, the prominent molecular absorption 
bands seen in the best-fit E60 (thin cloud) model are substan- 
tially reduced in the A60 (thick cloud) model, with AE60 lying 
in between. The planet’s flat 3-5 /im SED is best fit by the A60 
model. 

Although the / 2 minima for all four of the models we consider 
are sharply peaked, none yield fits falling within the 68% 
confidence interval. The fits from the E60 and AE60 models are 
particularly poor, ruled out at a greater than 5a level, whereas 
the A60 model quantitatively does better but still is ruled out 
as an acceptably fitting model (C.L. ~ 3.9a). The best-fit 
AMES-DUSTY model fits the SED even better than the A60 


model, with parameters of 7 dT = 1700 and log(g) = 3.5 and 
a radius of r — 1.35 Rj, similar parameters to those found in 
Bonnefoy et al. (2013). However, the best-fit DUSTY model 
still falls outside the 68% confidence limit (C.L. = 0.84). These 
exercises suggest that the atmospheric parameters assumed in 
the models need to be modified in order to better reproduce the 
/I Pic b photometry. To achieve this goal, we restrict ourselves 
to thick clouds and consider more carefully the impact of dust 
size. 

4.2. A4, Thick Cloud/ Small Dust Models 
4.2.1. The Effect of Small Dust Particles 

Our analyses in the previous section show the extreme mis- 
match between standard L dwarf atmosphere models assuming 
thin clouds and large dust particles and the data. While our y 2 
values for the Burrows thick-cloud, large-dust-particle models 
are systematically much lower, they likewise are a poor match 
to the data. In contrast, fits from the AMES-DUSTY models 
only narrowly lie outside the 68% confidence interval. 

A closer inspection of the best-fitting models in each case 
(right-hand panels of Figure 9) illustrates how they fail. The 
main difficulty with matching these models to the p Pic b 
spectrum is the planet’s flat SED from 2 /xm to 5 // m , where 
models tend to underpredict the flux density at 3.1 /xm and/or 
M' . The slope from J to K s is also a challenge. Reducing dust 
sizes can further fill in absorption troughs by increasing the 
opacities of the clouds. The AMES-DUSTY model, however, 
appears to overcorrect as its spectrum exhibits sharp peaks due 
to its submicron sized grains that degrade its fit to the data. 
Therefore, we consider grain sizes intermediate between those 
in A60 and AMES-DUSTY (i.e., -1-30 /xm). 

A4 thick-cloud, small-dust-particle models. As the primary 
difference between these models is the typical/modal particle 
size, we introduce here a new set of atmosphere models with 
the same A-type, thick cloud assumption but with modal particle 
sizes slightly larger than those characteristic of dust in the 
AMES-DUSTY models but significantly smaller than previous 
Burrows models. We nominally adopt 4 /xm as our new modal 
particle size, comparable in wavelength to the peak flux density 
of p Pic b in F v units. Thus, we denote these models as “A4,” 
thick-cloud, small-dust-particle models. 

Figure 10 illustrates the effect of dust on the planet spectrum 
for modal particle sizes of 3, 5, 30, and 50 /xm and a temperature 
and surface gravity consistent with those expected to reflect 
P Pic b based on planet cooling models (7’ c f f = 1600 K, 
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Figure 9. / 2 distributions (top left) and best-fitting models (top right, middle panels) for our limiting cases: the E60 thin cloud/large dust models, the AE60 
moderately thick cloud/large dust models, and the A60 thick cloud/large dust models. Bottom panels: the y 2 distributions and best-fit AMES-DUSTY models that 
assume ISM-sized dust grains. The horizontal lines in the x 2 panels display the 68% and 95% confidence limits. The pink bars roughly denote the model-predicted 
flux densities at the filters’ central wavelength positions. We identify the passbands along the bottom of the top-right panel. 

(A color version of this figure is available in the online journal.) 


l°g(g) = 3.8-4, r ~ 1.5 Rj', Burrows et al. 1997; Baraffe 
et al. 2003; Lagrange et al. 2010; Spiegel & Burrows 2012; 
Bonnefoy et al. 2013). As particle sizes decrease, the water 
absorption troughs at 1.8 /tm and 2.5 /im weaken. Likewise, the 
deep absorption trough at ~3.3 jtm and 4.5 // m that is usually 
diagnostic of carbon chemistry (e.g., Hubeny & Burrows 2007; 


Galicher et al. 2011) is filled in. Overall, the spectrum flattens 
and becomes redder (shorter wavelength emission originates at 
higher altitudes), with weaker emission and a steeper slope at J 
to K s . This reddening explains the difference in best-fit effective 
temperature between the AMES-DUSTY model and the 60 /im 
dust models. 


13 


The Astrophysical Journal, 776:15 (19pp), 2013 October 10 


Currie et al. 



Wovelength(Microns) 


Figure 10. Effect of atmospheric dust particle sizes on the shape of a massive 
planet SED. Here, the (a) = 50, 30, and 5 fim models depict log(g) = 4, 7j.fi = 
1600 K models while the (a) = 3 //ni model assumes log(g) = 3.8 and is scaled 
to match the luminosity of the (a) = 5 (im model. For small particle sizes, 
surface gravity signatures are weak and this parameter’s effect is primarily to 
change the flux scaling. 

(A color version of this figure is available in the online journal.) 

4.2.2. Model Fitting Procedure 

We follow the steps outlined in Currie et al. (2011a), where we 
perform two runs: one fixing the planet radius to the Burrows 
et al. (1997) hot-start predictions for a given T e g and log(g) 
and another where we consider a range of planet radii (as in 
the previous section). For the fixed-radii modeling, the 68% 
and 95% confidence limits now lie at / 2 = 7.01 and 12.6, 
respectively, whereas they are at 5.87 and 1 1.06 for the varying- 
radii fits as before. Similar to the Burrows A/E60 model runs, we 
consider a range of temperatures between 1400 K and 1900 K. 
To explore whether or not the fits are sensitive to surface gravity, 
we consider models with log(g) = 3.6, 3.8, 4, and 4.25. For the 
age of /3 Pic (formally, 8 to 20 Myr), this surface gravity range 
fully explores the masses (in the hot-start formalism) allowed 
given the radial velocity dynamical mass limits (Lagrange et al. 
2012b). 


To further explore the effect that carbon chemistry may have 
on our planet spectra, we take the best-fitting model from the 
above exercise, significantly enhance the methane abundances 
over solar and re-run a small grid of temperatures, to determine if 
departures from solar abundances yield a wider range of accept- 
able atmosphere parameters. Because variations in molecular 
abundances affect the depths of molecular absorption bands, we 
expect that such variations may improve our fit. 

4.2.3. Results 

Figures 1 1 and 1 2 and Table 8 present our results for fitting 
the Pic b data with the A4, thick cloud/small dust models. 
Quantitatively, these models better reproduce the Pic b SED. 
Fixing the planet radius to values assumed in the Burrows 
et al. (1997) planet cooling model, we find one atmosphere 
model — log(g) = 3.8, T e ff = 1600 K — consistent with the data 
to within the 68% confidence interval. A wide range of models 
are consistent with the data at the 95% confidence limit, covering 
± 0.2 dex in surface gravity and ±100 K in temperature. 

We can slightly improve upon these fits if we allow the 
planet radius to vary freely. In this case, the best-fitting models 
yield a slightly higher surface gravity of log(g) = 4-4.25 but 
the same temperature of 1600 K. But in contrast to the fixed- 
radius case above, a wide range of models are consistent with 
the data at the 68% confidence limit. In particular, all surface 
gravities considered in our model grid are consistent with the 
data provided that the temperature is 1600 K and the radius 
is rescaled accordingly: log(g) = 3.6-4.25, 7^ = 1600 K. 
Another set of models with a full range of surface gravities and 
a 250 K spread in temperature (1500-1750 K) is marginally 
consistent with the data. 

The methane-enhanced models are shown in Figure 13 for 
log(g) = 4 and T eS = 1575-1650 K. The 1575 K and 1600 K 
models (Figure 13) likewise produce good fits to the data (x 2 = 
5.13-5.3), where the 1650 K model barely misses the 68% 
cutoff. Thus, while best-fitting solar abundance models appear 
narrowly peaked at 7j, fT = 1600 K, the range in temperature 
enclosing the 68% confidence interval is larger when non-solar 
abundances are considered, at least extending from 1575 K to 
almost 1650 K. Changes in molecular abundances, as expected. 
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Figure 11. x 2 distribution (left) and best-fitting model (right) for the thick cloud, small dust A4 models, showing that we can achieve statistically significant fits to the 
data provided that the clouds are thick and the atmospheric dust particles are significantly smaller than those we previously assumed in matching L dwarf spectra (cf. 
Burrows et al. 2006). This model fitting ties the planet radius to predictions for hot-start models from Burrows et al. (1997). The horizontal dashed lines identify the 
95% confidence limit (top) and the 68% confidence limit (bottom). 

(A color version of this figure is available in the online journal.) 
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Figure 12. Same as Figure 11 except treating the planet radius as a free parameter. Here we show the log(g) = 4, T e ff = 1600 K model, although the synthetic 
spectrum’s appearance and its agreement with the data at neighboring gridpoints in surface gravity is nearly identical. 

(A color version of this figure is available in the online journal.) 
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Figure 13. SED fits adopting the nominally good-fitting atmosphere parameters depicted in Figure 12 but enhancing the atmosphere of methane. 

(A color version of this figure is available in the online journal.) 


allow us to very slightly improve the SED fit. However, thick 
clouds and small dust grains are likely still needed to match 
the emission from p Pic b, since given molecules (e.g., CH4) 
by themselves do not change fluxes comparably at 1-3 /xm and 
3-5 /xm. 

In summary, adopting the Burrows et al. (1997) hot-start 
models to set our planet radii and the A4 thick cloud/ small 
dust atmosphere models, we derive log(g) = 3.8 and 7/ f f = 
1600 K for p Pic b. Allowing the radius to vary and considering 
non-solar carbon abundances, we derive log(g) = 3.6-4.25 and 
T ef f = 1575-1650 K, meaning that the planet temperature is well 
constrained but the surface gravity is not. However, in Section 5 
we narrow the range of surface gravities to log(g) = 3.8 ± 0.2, 
as higher surface gravities imply planet masses ruled out by 
dynamical estimates. 

4.2.4. Varying Grain Sizes and Fits Over Other 
Model Parameter Spaces 

The models considered in the previous subsections assume 
thick clouds, dust grains with a modal size of 4 /xm, and (in 
most cases) solar abundances. Although we achieve statistically 
significant fits to the ft Pic b photometry with these models, 
our exploration of model parameter space is still limited. While 
an exhaustive parameter space search is beyond the scope of 
this paper, we argue here that models of either thick clouds or 
small dust grains are unlikely to be good-fitting models. Thus, 
small grains and thick clouds are likely important components 
of p Pic b’s atmosphere required in order to fit the planet’s 
spectrum. 


To consider the robustness of our results concerning the modal 
grain size, we also ran some model fits for modal particle sizes 
of 3, 5, 10, and 30 ptm. The models with 3 and 5 /xm modal sizes 
yielded fits slightly worse than those with modal sizes of 4 /xm. 
For example, models with modal sizes of (a) = 3 and 5 /xm, 
7/^ = 1600 K, log(g) = 3.8, and a freely varying planet radius 
yield x 2 of 6.31 and 6.28, respectively. These values lie slightly 
outside the 68% confidence interval, although they are still 
smaller than those from the best-fit DUSTY models. In contrast, 
models with (a) = 10 and 30 /xm fit the data significantly worse 
( x 2 = 10.0 and 19.6, respectively). 

Similarly, our investigations show that small dust grains do 
not obviate the need to assume thick, A-type clouds in our 
atmosphere models. For example, adopting the AE-type cloud 
prescription, modal particle sizes of 5 /xm, a temperature of 
T e ff = 1600 K, and a surface gravity of log(g) = 3.8-4, our 
model fits are substantially worse than the A4-type models and 
even the AMES-DUSTY models and are easily ruled out (x 2 ~ 
15-40). The AE-type cloud prescription fails to reproduce the 
P Pic b spectrum because by confining clouds to a thinner layer 
the r = 1 surface varies too much in and out of molecular 
absorption features such as CH4 and CO. In disagreement with 
the p Pic b SED, the AE model spectra thus have suppressed 
emission at «3 and 5 /xm and an overall shape looking less like 
a blackbody. 

In contrast, non-solar abundances may slightly widen the 
range of parameter space (in radius, temperature, gravity, etc.), 
yielding good fits. The methane-rich model from the previous 
section adopting (a) = 5 /xm instead of 4/xm, log(g) = 3.8, 
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Figure 14. Left: luminosity evolution for hot-start models from Baraffe et al. (2003), comparing ft Pic b’s luminosity as derived in this work with that for the directly 
imaged planets HR 8799 bcde, the planet/low-mass brown dwarf k And b, and other very low mass brown dwarf companions. The solid circle denotes ft Pic b’s 
nominal positions, whereas the open circle identifies its effective position on this plot if it formed after 5 Myr. Right: radius evolution for planets with masses of 
1-10 Mj for the “hot-start” and “cold-start” planet cooling models from Spiegel & Burrows (2012). The radius error bars define the range of planet radii from models 
consistent with the data to within the 68% confidence limit. 

(A color version of this figure is available in the online journal.) 


and T e g = 1600 K still yields a fit in agreement with the data 
to within the 68% confidence limit (/ 2 = 5.59). Thus, within 
our atmosphere modeling approach, we need (1) grains several 
microns in size, comparable to the typical sizes of grains in 
debris disks and (2) thick clouds to yield fits consistent with 
the data to within the 68% confidence limit. These results are 
not strongly sensitive to chemical abundances although varying 
the range of abundances may slightly widen the corresponding 
range of other parameter space (in temperature, gravity, etc.), 
yielding good-fitting models. 

5. PLANET RADII, LUMINOSITIES, MASSES, 

AND EVOLUTION 

From the set of models that reproduce the Pic b SED 
in the 68% confidence interval, we derive a range of planet 
radii, luminosities, and inferred masses. The planet radii for 
each model run are given in Table 8. Interestingly, all of our 
la solutions fall on or about R ~ 1.65 Rj with very little 
dispersion (±~0.05 dex). If we consider the range of radii for 
a given atmosphere model consistent with the data to within 
the 68% (or 95%) confidence interval regardless of whether the 
given radius is the best-fit one, then the range in acceptable 
radii marginally broadens: r = 1.65 ± 0.06. Note that these 
radii are larger than those inferred for HR 8799 bcde based 
either on their luminosities and hot-start cooling tracks (Marois 
et al. 2008, 2010b) or from atmosphere modeling, where in 
Currie et al. (201 la) and Madhusudhan et al. (2011) our best-fit 
models typically had R ~ 1.3 Rj. The range in inferred planet 
luminosities is even narrower. The values inferred from our 
best-fit models center on log (L/L Q ) = —3.80 with negligible 
intrinsic dispersion (± 0.01 dex). The uncertainty in ( J > Pic’s 
distance affects both our radius and luminosity determinations. 
Treating the distance uncertainty (± 1 pc) as a separate, additive 
source of error, /I Pic b’s range in radii is 1.65 ± 0.06 Rj and its 
luminosity is log (L/L Q ) = —3.80 ± 0.02. 

From our best-fit surface gravities and inferred radii, we can 
derive the mass of the planets inferred from our modeling. 
Adopting the hot-start formalism without rescaling the radius, 
our modeling implies a best-fit planet mass of ~7 Mj ; the range 
covering the 95% confidence limit in 5-9 Mj. If we allow the 
radius to freely vary, we derive a range of 4 Mj to 18.7 My, 
where the spread in mass reflects primarily the spread in surface 


gravity from best-fitting models (log(g) = 3.6-4.25). However, 
RV data limit (3 Pic b’s mass to be less than 1 5 Mj if its semimajor 
axis is less than 10 AU, which appears to be the case (Lagrange 
et al. 2012b; Chauvin et al. 2012; Bonnefoy et al. 2013). Thus, 
limiting the atmosphere models to those whose implied masses 
do not violate the RV upper limits (ones with log(g) = 3.6-4), 
our best-estimated (68% confidence limit) planet masses are 
~7* My. 

Planets cool and contract as a function of time and we can 
compare our inferred luminosities and radii with planet cooling 
models. Figure 14 compares the inferred planet luminosity to 
hot-start planet evolution models from Baraffe et al. (2003). For 
context, we also show the luminosities of other 5-100 Myr 
old companions with masses that (may) lie below 15 My: 
GSC 06214 B, 1RXJ 1609 B, HR 8799 bcde, AB Pic B, and 
k And b. From our revised luminosity estimate, the Baraffe et al. 
(2003) hot-start models imply a mass range of ~8 — 12 Mj if the 
planet’s age is the same as the star’s inferred age ( 1 2^ 8 4 Myr; 
Zuckerman et al. 2001). If we use the Burrows et al. (1997) 
hot-start models, we obtain nearly identical results of 9-13 Mj. 
These masses are slightly higher than most of the implied masses 
from our atmosphere modeling but still broadly consistent with 
them and with the dynamical mass upper limits of 15 Mj from 
Lagrange et al. (2012b). Note also that the luminosities and 
planet radii are completely inconsistent with predictions from 
low-entropy, cold-start models for planet evolution. 

Still, the right-hand panel of Figure 14 highlights one possible 
complication with our results, namely that our best-estimated 
planet radii are near the upper end of the predicted range for 
5-10 Mj companions in the hot-start formalism. For the hot- 
start models presented in Burrows et al. (1997) and Baraffe 
et al. (2003), 5-10 Mj companions are predicted to have radii 
of ~ 1 .5-1 .6 Rj. For the hot-start models presented in Spiegel & 
Burrows (2012), the predicted range for 5-10 Mj planets covers 
Ril.4-1.5M/. 13 

To reduce the planet radius of ~1.65 Rj by ~10% while 
yielding the same luminosity requires raising the effective 


13 This mismatch does not mean that the AMES-DUSTY models, whose fits 
to the data imply planet radii of ~ 1.3 Rj and lie just outside the 68% 
confidence limit, are preferable. The best-fit AMES -DUSTY radii lie below 
the radii predicted for 5-10 Mj objects at ft Pic b’s age and are only consistent 
for “warm-start” models that imply lower luminosities and colder temperatures 
than otherwise inferred from the AMES-DUSTY fits. 
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temperature from ~ 1 600 K to ~1700 K. This change is small 
and atmospheric modeling of p Pic b and similar substellar 
objects is still in its early stages. Thus, it is quite plausible that 
future modeling efforts, leveraging on additional observations 
of p Pic b and those of other planets with comparable ages and 
luminosities, will find quantitatively better fitting solutions that 
imply smaller planet radii and higher temperatures. We consider 
this statement to be the most likely explanation for our larger 
radii. 

Alternatively, we can bring the atmosphere modeling-inferred 
radius into more comfortable agreement with hot-start evolu- 
tionary models if p Pic b is ^7 Myr old or less. For a system 
age of »12 Myr, this age is consistent with the planet forming 
late in the evolution of the protoplanetary disk that initially sur- 
rounded the primary. Even adopting the lower limit on p Pic’s 
age (8 Myr), p Pic b may still need to be younger than the star. 
While most signatures of protoplanetary disks around 1-2 M Q 
stars disappear within 3-5 Myr, some ~ 1 0% -20% of such 
stars retain their disks through 5 Myr (Currie et al. 2009; Currie 
& Sicilia-Aguilar 2011; Fedele et al. 2010). Several 1-2 M e 
members of Sco-Cen and h and / Persei apparently have even 
retained their disks for more than 10 Myr (Pecaut et al. 2012; 
Bitner et al. 2010; Currie et al. 2007), comparable with or greater 
than the age of p Pic. Models for even rapid planet formation by 
core accretion predict that several Myr elapse before the cores 
are massive enough to undergo runaway gas accretion at p Pic- 
like separations (Kenyon & Bromley 2009; Bromley & Kenyon 
2011). 

In Figure 14, the open circles depict a case where p Pic b 
formed after 5 Myr, effectively making the planet 5 Myr younger 
than the star, where the implied masses and radii overlap better 
with our atmospheric modeling-inferred values. The overlap 
is even better for some hot-start models such as COND, which 
predict larger planet radii at s»5-10 Myr than depicted here. Note 
that a young p Pic b as depicted in Figure 14 with an implied 
mass of AO 5 Mj is still consistent with a scenario where the 
planet produces the warped secondary disk (cf. Dawson et al. 
2011). 

6. DISCUSSION 
6.1. Summary of Results 

This paper presents and analyzes new/archival VLT/NaCo 
and Gemini/NICI 1-5 p . m photometry for p Pic b. These data 
allow a detailed comparison between p Pic b’s SED and that of 
field brown dwarfs and other low-mass substellar objects such as 
directly imaged planets/candidates around HR 8799 and k And. 
Using a range of planet atmosphere models, we then constrain 
P Pic b’s temperature, surface gravity and cloud properties. Our 
study yields the following primary results. 

1. The near-IR ( JHK S ) colors of p Pic b appear fairly con- 
sistent with the held L/T dwarf sequence. Compared with 
other young, low-mass substellar objects, p Pic b’s near- 
IR colors bear the most resemblance to late M to early T 
dwarfs such as Luhman 16B and k And b. From its near-IR 
colors and color-magnitude positions, p Pic b’s near-IR 
properties most directly resemble those of a L2-L5 dwarf. 

2. p Pic b’s mid-IR properties imply a significant departure 
from the held L/T dwarf sequence. The planet is slightly 
overluminous at L' and signihcantly overluminous at M', 
with deviations from the held L dwarf sequence matched 
only by GSC 062 14B and k And b. The mid-IR portion of 
P Pic b’s SED appears more like that of a late L dwarf or 


low surface gravity mid L dwarf. The broadband JHK S L' 
photometry for p Pic b also closely resembles that of k 
And b. However, it is unclear whether any object matches 
P Pic b’s SED at all wavelengths for which we have 
measurements. Its 3. 1 /am brightness and 3.8-5 /im spectral 
shape are particularly difficult to match. 

3. Compared with limiting-case atmosphere models E60 
(large dust conhned to very thin clouds), AE60/A60 
(large dust conhned to moderately thick/thick clouds), 
and DUSTY (copious small dust everywhere in the atmo- 
sphere), p Pic b appears to have evidence for thick clouds 
consistent with a high T e ff and low surface gravity. We fail to 
hnd any E60/AE60/A60 model providing statistically sig- 
nificant hts over a surface gravity range of log(g) = 4-4.5 
and any 7/n. The DUSTY models come much closer to 
yielding statistically signihcant hts but mismatch the planet 
flux at J, K s , [3.1], and M' . From these hducial compar- 
isons, we infer that p Pic b’s atmosphere shows evidence 
for clouds much thicker than those assumed in the E60 
models but is slightly less dusty than the DUSTY models 
imply. 

4. Using thick cloud models with particle sizes slightly larger 
than those found in the interstellar medium (ISM; (a) = 
4 /im), we can match p Pic b’s SED in both the near- and 
mid-IR. Assuming planet radii appropriate for the Burrows 
et al. (1997) “hot-start” models, we derive log(g) = 3.80 
and 7/ if = 1600 K for p Pic b. Allowing the radius to freely 
vary leaves the surface gravity essentially unconstrained, 
where models consistent with the data at the 68% confi- 
dence limit include log(g) = 3.6 — 4.25 and T e ff = 1600 K. 
Considering departures from solar abundances and elimi- 
nating models that imply masses ruled out by dynamical 
estimates, the acceptably fitting range of atmosphere pa- 
rameters cover log(g) = 3.6-4 and 7’ c rr = 1575-1650 K. 

5. Using our best-fit atmosphere models and eliminating 
models inconsistent with p Pic b’s dynamical mass upper 
limit, we derive a mass of 7 Mj for a fixed radius and 
1% Mj for a scaled radius within the hot-start formalism. 
Our best- fit planet radius is ~ 1 .65 ± 0.06 Rj and our best-fit 
luminosity is log(L/L G ) = —3.80 ± 0.02. 

6. While our derived luminosity and radius for p Pic b rule 
out cold-start models, the radius is near the upper end of 
predicted radii for hot-start-formed planets of P Pic’s age. 
As the planet only needs to be ~ 1 00 K hotter to easily 
eliminate this discrepancy, this mismatch likely identifies 
a limitation of the atmosphere models. Alternatively, if 
P Pic b has a significantly younger age than the star’s age, 
consistent with it forming late in the protoplanetary disk 
stage, our derived radius is comfortably within the range 
predicted by hot start models. 

6.2. Comparisons with Other Recent p Pic b Studies 
6.2.1. Currie et al. 2011b 

In our first-look analysis of the atmosphere of p Pic b (Currie 
et al. 201 lb), we compared its K„, L', and [4.05] photometry with 
an array of atmosphere models, from atmospheres completely 
lacking clouds to those with the Model A-type thick clouds that 
extend to the visible surface of the atmosphere. In that paper, 
we found that the AE thick cloud models from Madhusudhan 
et al. (201 1) yielded the smallest / 2 value. The hts degraded at 
about the same level for the Model A thick cloud and Model E 
“normal” L dwarf atmosphere prescriptions, while the cloudless 
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case fared the worst. Currie et al. (2011b) conclude that while 
the AE thick cloud model quantitatively produced the best fit, 
the existing data were too poor to say whether the clouds on 
P Pic b were any different in physical extent, in mean dust 
particle size, etc. from those for field L dwarfs with the same 
range of temperatures. 

Our present study greatly improves upon the analyses in 
Currie et al. (2011b). First, our photometry covers seven pass- 
bands, not three, over 1.25-4.8 /im, not 2.18-4.05 /xm. This ex- 
panded coverage allows far firmer constraints on p Pic b’s atmo- 
spheric properties. In particular, our new photometry strongly 
favors the Model A thick-cloud prescription over the AE 
model, largely due to the relatively low planet flux densities 
at 1.25-1.65 /x m and the relatively high flux densities at 3.1 /im 
and M' / 4.8 /im, trends that the Model A cases consistently re- 
produce better. While all models considered in Currie et al. 
(2011b) assumed a modal particle size of 60 /im for dust en- 
trained in clouds, our fits improve if we use smaller dust sizes. 
The combined effect of thicker clouds and smaller particle sizes 
favors atmosphere models with a slightly higher surface gravity 
and temperature than the best-fit model in Currie et al. (201 lb). 
Our new data more clearly demonstrate the failure of the E 
models successful in fitting most of the field L dwarf sequence 
and thus better distinguish p Pic b’s atmosphere from that of a 
typical cloudy field L dwarf. 

6.2.2. Bonnefoy et al. (2013) 

Bonnefoy et al. (2013) presented new photometry for p Pic b 
in the J , H, and M' filters from data taken in 2011 and 2012. 
The J and H detections are firsts and greatly expand the wave- 
length coverage for p Pic b’s SED. Their M' detection is the 
first well-calibrated detection, building upon and following the 
detection presented in Currie et al. (2011b), which lacked con- 
temporaneous flux-calibration data to provide precise photom- 
etry. Bonnefoy et al. (2013) then combined these measurements 
with their previously published K s and L' photometry and [4.05] 
from Quanz et al. (2010). 

In general, our study clarifies and modifies, instead of 
contradicting, the picture of p Pic b constructed in Bonnefoy 
et al. (2013). In the same data sets, the S/N of our detections is 
slightly higher but our photometry agrees within theirs derived 
from their CADI, RADI, and LOCI reductions within their 
adopted photometric uncertainties (~0.2-0.3 mag). We derive 
smaller photometric uncertainties, owing to a more uniform 
throughput as a function of azimuthal angle, probably due to 
our pixel masking technique and SVD cutoff in A-LOCI (see 
also Marois et al. 2010a). We concur that the planet’s mid-IR 
colors are unusually red and highlight a potentially strong, new 
disagreement with field L dwarfs at 3.1 /xm. 

We agree with Bonnefoy et al.’s general result that the best- 
fitting atmosphere models are those intermediate between the 
AMES-DUSTY models (submicron-sized dust everywhere) and 
the COND or BT-Settl models (no dust/clouds or very thin 
clouds). Quantitatively, the y 2 values we derive are much larger 
than the best-fitting models in Bonnefoy et al. (i.e., 7 vs. 3 
for AMES-DUSTY) because our photometric uncertainties are 
significantly smaller. Our analyses point to thick clouds and 
particle sizes small compared with the range typically used in 
the Burrows et al. (2006) models but larger than the ISM-like 
grains in the AMES-DUSTY models. The temperatures, surface 
gravities, and luminosities they derive are generally consistent 
with our best-fit values. 


While Bonnefoy et al. (2013) derive a lower limit to the initial 
entropy of 9.3 kB baryon 1 , we do not provide a detailed similar 
analysis since the inferred entropy range depends on the planet 
radius which, considering our studies together, is very model 
dependent. Similarly, entropy depends on the planet mass (for 
which there still is some range) and the planet’s age (which is 
very poorly constrained). Still, we agree that cold start models 
are ruled out for p Pic b as they fail to reproduce the inferred 
luminosity and radii of the planet determined from both our 
studies. 

6.3. Future Work to Constrain p Pic b’s Properties 

Deriving p Pic b’s mass and other properties is difficult since 
these parameters are based on highly uncertain parameters such 
as the planet’s age and its entropy at formation. However, 
dynamical mass limits can be derived from continued RV 
measurements (Lagrange et al. 2012b). As these limits depend 
on p Pic b’s orbital parameters, future planet astrometry may be 
particularly important in constraining p Pic b’s mass. If p Pic b 
is responsible for the warp observed in the secondary debris disk 
(Golimowski et al. 2006), planet-disk interaction modeling can 
likewise yield a dynamical mass estimate (Lagrange et al. 2009; 
Dawson et al. 2011), provided the planet’s orbit is known. 

Finally, while our models nominally assume solar abun- 
dances, we showed that changing the methane abundance might 
yield marginally better fits to the data. Near-infrared spectro- 
scopic observations of p Pic b, as can be done soon with GPI 
and SPHERE, may clarify its atmospheric chemistry. Future 
observations with GMTNIRS on the Giant Magellan Telescope 
should be capable of resolving molecular lines in p Pic b’s at- 
mosphere (Jaffe et al. 2006), providing a more detailed look at 
its chemistry, perhaps even constraining its carbon to oxygen 
ratio and formation history (e.g., Oberg et al. 2011; Konopacky 
et al. 2013). 
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